Helicity conservation in nonlinear mean-field solar dynamo 
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We explore the impact of magnetic helicity conservation on the mean-field solar dynamo using the 
axisymmetric dynamo model which includes the subsurface shear. Our results support the recent 
findings by Hubbard & Brandenburg [6], who suggested that the catastrophic quenching in the 
mean-field dynamo is alleviated if conservation of the total magnetic helicity is taken into account. 
We show that the solar dynamo can operate in the wide rage of the magnetic Reynolds number up to 
10 6 . We also found that the boundary conditions for the magnetic helicity influence the distribution 
of the a-effect near the solar surface. 

Keywords: Turbulence: Mean-field magnetohydrodynamics; Sun: magnetic field; Stars: activity 
- Dynamo 

PACS numbers: 



INTRODUCTION 



The basic idea for the solar dynamo action was de- 
veloped by Parker [13]. He suggested that the toroidal 
component of the magnetic field of the Sun is stretched 
from the poloidal component by the differential rotation 
(f2 effect) and the cyclonic motions (a effect) return the 
part of the toroidal magnetic field energy back to the 
poloidal component. This is the so-called afi scenario. 
This mechanism is implemented in the wide range of the 
solar dynamo models (see review bySJ- 

The effect of the turbulence in the mean-field dynamo 
is represented by the mean electromotive force £=u x b, 
where u and b are the fluctuating velocity and mag- 
netic fields. In the simplest case it can be found that 

£ = a B + P ' x B^ — t\t (V x B) , where ao is the a 

— (p) 

effect, V is the turbulent pumping and r\x is the turbu- 
lent diffusivity [TU] . The a effect is a pseudo-scalar (lacks 
the mirror symmetry) which is related to the kinetic he- 

T c 

licity of the small-scale flows, i.e., ao = tt u ' ^ Xu : 

o 

where r c is correlation time of turbulent motion. Pou- 
quet et al. [19] showed that the a effect is produced not 
only by kinetic helicity but also by current helicity, and 



it is ao = — - 



u • Vxu - 



b - V x b 



The latter ef- 
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feet can be interpreted as resistance of magnetic fields 
against to the twist by helical motions. It leads to the 
concept of the catastrophic quenching of the a effect by 
the generated large-scale magnetic field. It was found 
a (0) 



that ao (B) 



j, where R m is magnetic 



1 + R m [B/ B eq ) 
Reynolds number (see, Eland references therein). In case 
of R m ^> 1, the a effect is quickly saturated for the 
large-scale magnetic field strength that is much below the 



equipartition value B eq ~ ypfi u 2 . The result was con- 
firmed by the direct numerical simulations (DNS) |12j . 
The catastrophic quenching (CQ) is related to the dy- 
namical quenching of the a effect. It is based on conser- 
vation of the magnetic helicity, \ = ab (a is fluctuating 
part of the vector potential) and the relation between the 
current and magnetic helicities he = b • V X b ~ x/^ 2 , 
which is valid for the isotropic turbulence |11|. The evolu- 
tion equation for % can be obtained from equations that 
govern a and b, it reads as follows [3]: 

% = -2(£.B)-^--V-^-»jS-J, (1) 

at R m T c 

where, in following to Kleeorin & Rogachevskii (jJJ, we in- 
troduce the helicity fluxes T x = axuxB — ax (uxb). 
The helicity fluxes are capable to alleviate the catas- 
trophic quenching (see, e.g.. [31 12^1 and references therein). 
Generally, it was found that the diffusive fluxes, which 
are ~ 7y x Vx, where r\ x is the turbulent diffusivity of the 
magnetic helicity, work robustly in the mean-field dy- 
namo models but it requires r\ x > t\t to reach \B\ > 
0.1B eg . 

Another possibility to alleviate the catastrophic 
quenching is related with the non-local formulation of 
the mean-electromotive force [2]. In fact, the Babcock- 
Leighton (BL) type dynamo is the the special case of the 
mean-field dynamo with the nonlocal a effect. Kitchati- 
nov & Olemskoy [8] found that the nonlocal a effect and 
the diamagnetic pumping can alleviate the catastrophic 
quenching. The results by Brandenburg & Kapyla [1 
show that the strength of the quenching can depend on 
the model design. Therefore, the problem of the catas- 
trophic quenching is actual for different types of the solar 
dynamo. 

Recently, Hubbard & Brandenburg 6 j revisited the 
concept of catastrophic quenching and showed that for 
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the shearing dynamos the Eq. ([!]) produces the nonphysi- 
cal fluxes of the magnetic helicity over the spatial scales. 
They suggested to cure the problem using the global con- 
servation law for the total magnetic helicity that can be 
written as follows: 



dt 



{x + A-B}dV = - v / {B-J + b-i}dV(2) 



V-T x dV 



where integration is done over the volume that comprises 
the ensemble of the small-scale fields. We assume that 
T x is the diffusive flux of the total helicity which is re- 
sulted from the turbulent motions. Ignoring the effect of 
the meridional circulation we write the local version of 
the Eq.Q as follows [5]: 



dx _ d(A-B) 
dt ~ 



dt 



ryB ■ J V-F x 



(3) 



where J- x = — r\ x V (x + A • B) . In the paper we employ 
this equation for the dynamical quenching of the a effect 
in the solar dynamo model and show how it works in 
the range of the magnetic Reynolds number R m = 10 3-6 
those are typical for the astrophysical conditions. 



BASIC EQUATIONS 

We study the mean-field induction equation in a per- 
fectly conducting medium: 



dB 

~dt 



V x (S + U x B) , 



(4) 



where £ = u X b is the mean electromotive force, with 
u, b being fluctuating velocity and magnetic field, re- 
spectively, U is the mean velocity (differential rotation), 
and the axisymmetric magnetic field is: 



B = e^B + V x 



Ae,, 



where 9 is the polar angle. The mean electromotive force 
£ is given by Pipin [TJ]. It is expressed as follows: 



£i = 



■nr.. 



(A) 



Bj-riijkVjBk. 



The tensor ctij represents the a-effect, 7y is the turbu- 
lent pumping, and r]ijk is the diffusivity tensor. The a 
effect includes hydrodynamic and magnetic helicity con- 
tributions, 



(5) 



n ■ 2 a (H) . (M) 
atij = C a sin 9alj + 



(6) 



The details in expressions for the kinetic part of the a 

(H) „ (A 

*V , as well as 7^ 



effect ot[j , as well as 7^ and rjijk can be found in 



The contribution of small-scale magnetic helicity x = 
a • b (a is the fluctuating vector-potential of the magnetic 
field) to the a-effect is defined as 



(M) _ of (a) 



XTc 
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XTc 



(7) 



The nonlinear feedback of the large-scale magnetic field 
to the a-effect is described by a dynamical quenching due 
to the constraint of magnetic helicity conservation given 
by Eq.(|3]). The effect of turbulent diffusivity, which is 
anisotropic due to the Coriolis force, is given by: 

Vijk = 3??t {f2/i W - / 2 (d) ) e ijk - 2/ 1 (a) e i e„e njfc (8) 



+Csl 



where functions /j"'^ 4} depend on the Coriolis number. 
They can be found in Pipin [T5] . The last part of Eq. ^ 
stands for the Q x J effect 20 . The DNS dynamo ex- 
periments support the existence of the dynamo effects 
induced by the large-scale current and global rotation 

We matched the potential field outside and the perfect 
conductivity at the bottom boundary with the standard 
boundary conditions. For the magnetic helicity we em- 
ploy x — at the bottom of the convection zone. At the 
top we use two types of the boundary conditions like 



v r x = 0, 

V r (x + A-B) = 0. 



(9) 
(10) 



To evolve the Eq.([3| we have to define the large-scale 
vector potential for each time-step. For the axisymmet- 
ric large-scale magnetic fields where the vector-potential 
is A = e^A/ (r sin 9)+re r T. The toroidal part of the vec- 
tor potential is governed by the dynamo equations. The 
poloidal part of the vector potential can be restored from 
equation V x (rT) = e^B. The restoration procedure is 
trivial for the pseudo-spectral numerical schemes which 
are based on the Legendre polynomial decomposition for 
the latitude profile of the large-scale toroidal field. The 
choice of parameters in the dynamo is justified by our 
previous studies [15]. Here we use C a = 0.03, Cs = \C a 
and the diffusivity dilution factor C n = 0.1. The pa- 
rameters of the models are summarized in the Table [U 

We estimate the turbulent parameters in the model on 
the base of the mixing-length theory and use as the ref- 
erence the solar interior model computed by Stix |23j . 
The differential rotation profile is like that suggested by 
Pipin & Kosovichev [17] . Fig. [I] shows the radial pro- 
files of the a effect components and profiles of the back- 
ground turbulent diffusivity C v t]t, the isotropic, rf 1 ^, 
and anisotropic, r)^ A ', parts of the magnetic diffusivity as 
well as profile for the il x J effect. To quantify the mirror 
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Figure 1: Left, (a), the profiles of the a effects components for 
the 8 = 45°. Right, (b), the profiles of the background tur- 
bulent diffusivity C v r]r, the isotropic, rf 1 ', and anisotropic, 
r/ A ', parts of the magnetic diffusivity and fl X J effect (also 
known as 8 effect |20| . 



Figure 2: Left panel (a) shows the ratio of the maxima 
of the toroidal magnetic field strength and the equipartition 
value with time; (b) shows the parity indexes evolution in 
the models. We employ the running average to filter out the 
separate cycles. 



Table I: Summary of parameters of the models. Here, QT 
stands for the quenching type and BC - the boundary condi- 
tions. 

Model QTli, QT2i QT1 2 , 3 QT2 2 , 3 



R* 



10° 



QT Eq. (131 ,Eq. < 3 1 Eq. (31 



Vx 
BC 



10 "rjr 
Eq. |9l 



10 a ,10 b 



10 d ,10 s 
^qTfl 



Eq. H Eq. |9l,Eq. (lOl 



symmetry type of the toroidal magnetic field distribution 
relative to equator we introduce the parity index P: 



P = 



Eq — Ed 
E„ + Ed 



E d = j {B(r ,9)-B(r Q ,ir-9)) 2 sm8d9, 

E q = J (B (r , 9) + B (r , vr - 9)f sin 8d0, 

where Ed and E q are the energies of the dipole-like and 
quadruple-like modes, r — O.9i? . We define the dy- 
namical quenching type 1 (QT1) to be govern by Eq.(l I, 
and the dynamical quenching type 2 (QT2) - by Eq. (3 1. 

RESULTS 



Fig. [2] shows the long-term evolution of the maximum 
of the large-scale magnetic field strength in the convec- 
tion zone and the parity in the models. The energy of 
the toroidal magnetic fields in all models show the expo- 
nential grow in the beginning phase, which has duration 
about 10 diffusive time of the system. The greater ini- 
tial magnetic field strength, the shorter duration of the 
exponential grow phase. Two models QTli and QT2i 
have rather small diffusive fluxes of the helicity and the 
high magnetic Reynolds number R m = 10 6 . We consider 
them as the references. It is shown below, see Figures 
[2] and [4] that the model QT2i is not subjected to the 
catastrophic quenching while the model QTli saturates 
the toroidal magnetic field strength to the level that is 



much below the equipartition. The moderate diffusive 
flux with r\ x — 0.01?7t (model QTI3) does not alleviate 
the problem if the magnetic helicity evolution is governed 
by the Eq. Q. In the models QT1 2 and QT2 li3 the sat- 
uration level of the toroidal magnetic field strength is 
about 0.3i? O q, where B eq is the equipartition level of the 
magnetic field strength. It is about 0.5_B eg in the model 
QT2 2 which has R m = 10 3 . The saturation level in the 
QT2 types solar dynamo can be higher for the the greater 
C a . This question needs a separate study. In the model 
QTli as well as in all the QT2 models the parity index 
saturates to the dipole-like solution. In the QTl2,3 mod- 
els the asymptotic stage is not clear and they need much 
longer runs. 

The origin of difference in behaviour of the magnetic 
helicity evolution in the models with QT1 and QT2 has 
been discussed recently by Hubbard & Brandenburg @. 
Taking into the dynamo equation Q, the corresponding 
equation for the large-scale vector potential: 



OA 

~dt 



S + U x B, 



(11) 



where we assume the Coulomb gauge, we find the equa- 
tion which governs the large-scale helicity evolution: 



d (A • B) 

dt 



= 2S ■ B+V • ((fx A)-A x (U x B))(12) 



Therefore, Eq. ^ can be rewritten in form of Eq.Q: 



dx 
dt 



X 



-2 (£ ■ B) 

1 tj m 1 c 

-rjB ■ J - V • (£ x A 



V • ( Vx Vx) (13) 



The term (£ x A) consists of the counterparts of the 
sources magnetic helicity, which are represented by — IS ■ 
B, and the fluxes which result from pumping of the large- 
scale magnetic fields. The sources magnetic helicit y in 
the term —2 (£ ■ B) are partly compensated in Eq(13) 
by the counterparts in (5 x A) . This results in the spa- 
tially homogeneous quenching of the large-scale magnetic 



Figure 3: Snapshots of the magnetic field and helicity evolution inside the North segment of the convection zone for the model 
QT23. Panel (a) shows the field lines of the poloidal component of the mean magnetic field and the toroidal magnetic field 
(varies ±700G) by gray scale density plot. The bottom panel, (b), shows the large-scale (density plot) and small-scale magnetic 
helicity (contours) distributions. Both kinds of the magnetic helicity vary in the same range of magnitude. 



generation and alleviation of the catastrophic quenching 



problem. The last term in Eq( 12 ) contains the transport 



of the large-scale magnetic helicity by the large-scale flow. 

Fig. [3] shows the snapshots of the magnetic field and 
magnetic helicity (large- and small-scale) evolution in the 
North segment of the solar convection zone for the model 
QT23. We observe the drift of the dynamo waves which 
are related to the large-scale toroidal and poloidal fields 
towards the equator and towards the pole, respectively. 
The distributions of the large- and small-scale magnetic 
helicities show the correspondence in sign: positive to 
negative, and the other way around, respectively. This 
is in agreement with Eq. p]). It is seen that the neg- 
ative sign of the magnetic helicity follows the dynamo 
wave of the toroidal magnetic field. This can be re- 
lated to the so-called "current helicity hemispheric sign 
rule" (negative/positive sign of helicity dominate in the 
North/South hemisphere) which is suggested by the ob- 
servations (see Seehafer [55], Zhang et al. [55] and ref- 
erences therein). The origin of the helicity sign rule 
has been extensively studied in the dynamo theory (e.g., 
l5l [25]) . The similar snapshots for the QT1 model can be 
found in [16]. The main difference is that for the QT1 
model all the changes in the magnetic helicity evolution 
does not exactly follow to the dynamo wave inside the 
convection zone as it is demonstrated for the QT2 model. 

Fig. [4] shows variations of the radial profiles of the a 
effect and magnetic helicity with the cycle and the time- 
latitude diagrams for the dynamo model QT2 3 . For all 
the models the changes in the a effect arc concentrated 
to the top of the convection zone. This is due to the 
factor (pi 2 ) in relation between the current and mag- 
netic helicities. The difference in the boundary at the 
top results in different distributions of the a effect in the 
models QT2r and QT2 3 . We would like to notice that 
the QTI2 and QT2i have the same boundary conditions, 
however, the resulted distributions of the a effect differ 



drastically. The time-latitude diagrams for the dynamo 
model QT2 3 which are shown in Figure [4] are in qualita- 
tive agreement with observations. The same patterns are 
obtained in the models QT2 12 . We show the dynamical 
a effect as well. The model shows that with the bound- 



ary conditions given by Eq. (10 1 the a effect increases 



and has positive maxima at the growing phase of the cy- 
cle and it decreases, having the negative minima at the 
decaying phase of the cycle. 



CONCLUSION 

In the paper we studied the effect of the magnetic he- 
licity conservation in the mean-field solar dynamo model 
which is shaped by the subsurface shear. The results 
show that the solar dynamo can operate in the wide rage 
of the magnetic Reynolds number up to 10 6 if conserva- 
tion of the total magnetic helicity is taken into account. 

It was found that the boundary conditions for the mag- 
netic helicity influence the distribution of the a effect 
near the solar surface. For example, the dynamo wave 
becomes closer to equator when the diffusive flux of the 
total helicity is zero at the top of the convection zone 
because the dynamical increase of the a effect follows 
the dynamo wave (see Fig. |3j [4]). The situation is op- 
posite for the case when the diffusive flux is dominated 
by the large-scale helicity. Thus, the parity breaking in 
the solar dynamo can occur due to perturbations in the 
external boundary for the magnetic helicity. However, 
these points remain an open field for the future work. 
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Figure 4: First row, (a,b and c), shows variations of the a effect (a^ component) profiles with the cycle at colatitude 9 = 45°, 
for the model QT2i, QT23 and QTI2 (from left to right). The bottom, panel (d) shows the time - latitude variations of the 
toroidal field near the surface, r — 0.95-Rq, (contours ±300G) and the radial magnetic field at the surface (density plot) for 
the model QT23; panel (e) shows the same for the the toroidal field and the a effect (a^ component) (density plot). 
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